Time-reversal in dynamically-tuned zero-gap periodic systems 
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We show that short pulses propagating in zero-gap periodic systems can be reversed with 100% 
efficiency by using weak non-adiabatic tuning of the wave velocity at time-scales that can be much 
slower than the period. Unlike previous schemes, we demonstrate reversal of broadband (few cycle) 
pulses with simple structures. Our scheme may thus open the way to time-reversal in a variety of 
systems for which it was not accessible before. 
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Wave propagation in periodic systems has been the fo- 
cus of countless studies in several branches of physics. 
The most prominent of these systems are crystalline 
solids and their artificial analogues, known as photonic 
crystals (PhCs) phononic crystals etc.. The study 
of such systems is almost exclusively associated with 
the forbidden gaps, which are energy /frequency regimes 
where waves cannot propagate. These gaps are crucial 
to understanding the electronic and optical properties of 
solid-state systems; their artificial analogues have a vari- 
ety of possible applications, ranging from low-threshold 
lasers, light/sound guiding, filtering and switching [l| to 
medical ultrasound and nondestructive testing [H . 

Although the occurrence of bandgaps is generic to pe- 
riodic systems, in some special cases, the gaps can have a 
zero width. Similarly, some materials and structures sup- 
port two bands that cross symmetrically, without forming 
a gap. One of the most notable of those is graphene Q, 
or its photonic crystal analogue pj. Others are chiral 
metamaterials [3], transmission- line systems Q, bi-axial 
crystals @ and spin systems Q. Previous studies of 
such systems have already shown some unusual proper- 
ties [1,13, SB- However, these systems have been rela- 
tively unexplored, especially for pulse manipulations. For 
this purpose, the zero-gap system differs from a finite- 
gap system in the sense that while a pulse incident on a 
finite-gap system will mostly be reflected, a pulse incident 
on a zero-gap system will be almost perfectly admitted. 
Moreover, the proximity of the crossing bands allows for 
efficient transfer of energy between the bands using weak 
and slow modulations of the wave velocity. This enables 
light slowing-down, stopping or reversing etc.. 

In this Letter, we focus on time-reversal. A time- 
reversed pulse evolves as if time runs backward, thus 
eliminating any distortions or scattering that occurred 
at earlier times. It has applications in diverse fields such 
as medical ultrasound [9j, optical communications and 
adaptiv e op tics [Iol - fl2l ]. superlensing [l3j], ultrafast plas- 
monics [lj], biological imaging |l5j . THz imaging [l6| 
and quantum information and computing fl7j . In par- 
ticular, we study time-reversal of electromagnetic pulses 
which, to date, has been demonstrated for pulses of 
a relatively narrow spectrum or requires complicated 
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schemes [lOt I18H23I ] . Some of these techniques [21|,|22| are 
based on dynamically tuned PhCs. This novel subject of 
PhC research is experiencing a fast growth with appli- 
cations such as frequency shifting, switching and many 
others studied theoretically and experimentally [HI Ha l- 
In this Letter, we show that dynamically tuned zero-gap 
PhCs enable an effective and broadband reversal in con- 
ceptually simple systems which are well within contem- 
porary fabrication capabilities. Most importantly, unlike 
previous schemes, the zero-gap systems allow for reversal 
of pulses of unprecedented broad spectrum, thus, open- 
ing the way to reversal of few cycle pulses. In addition, as 
our scheme does not rely on any concept which is unique 
to optics, it can be applied for reversal of other wave sys- 
tems for which time-reversal was not accessible before. 

We consider the simplest example of a perfectly sym- 
metric zero-gap system, namely, a ID layered PhC satis- 
fying the quarter- wave stack (QWS) condition, i.e., when 
the indices and thicknesses of the two layers satisfy 0, [2(| 

nidi = n 2 d 2 . (1) 

Such a system has a zero-width gap at, e.g., 

A c = 2n\d\ = 2n 2 d 2 . (2) 

The band structure of such PhCs can be computed an- 
alytically [26|, see e.g. Fig. [TJa). At the crossing point 
(at zero Bloch momentum, K = 0), light can travel ei- 
ther in the forward or backward direction (i.e., on the 
positive/negative group- velocity band). An index mod- 
ulation causes a transfer of energy between these bands. 
Intuitively, accurate reversal can be achieved only if two 
conditions arc fulfilled. First, the modulation should be 
periodic in order to avoid any wavevector mixing, re- 
sulting in a vertical frequency-conversion, see Fig. [TJb). 
Second, the modulation should be much faster than the 
pulse duration or equivalently, the spectral content of 
the modulation should be much wider than that of the 
pulse; only such a non-adiabatic modulation ensures a 
frequency-independent frequency conversion, hence, an 
accurate reversal. Nevertheless, the modulation can still 
be much longer than the period of the pulse oscillations. 

In what follows, we analyze the QWS system and con- 
firm the intuitive arguments given above. However, we 
emphasize that the ideas and techniques we use here can 
be employed also in other zero-gap systems. 

Consider an electromagnetic pulsed plane-wave nor- 
mally incident on a QWS PhC (see Fig. [5J which is time- 
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FIG. 1: (Color online) (a) Reduced Brillouin zone plot of the 
band structure of a QWS PhC. The bands support pulse prop- 
agation in opposite directions (indicated by the red arrows), 
(b) Schematics of the frequency conversion process. 



modulated in the following form 

n(x,t) = nQwsix) + Ap(x)m(t), max[m(t)] = 1. 

Here, the static QWS riQy/s(x) and spatial pattern of 
the modulation p{x) have period d, and A is a constant 
representing the magnitude of the modulation. For a 
linear polarization, the Maxwell equations are given by 

E x (x, t) = -ti H t {x, t), H x = -e [n 2 (x, t)E(x, t)] ( , (3) 

or, if the 7J-field is eliminated, by 

E xx (x,t) = [v- 2 (x,t)E(x,t)] tt , (4) 

which is a one-dimensional scalar wave equation with a 
space and time-dependent velocity v = c/n(x,t). 
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FIG. 2: (Color online) Geometry of the pulse propagation 
through the QWS PhC. 

We proceed with an analysis similar to that first per- 
formed in [2^ . We define 

E(x,t) = n~^(x,t)[F(x,t)+B(x,t)], 
Z H(x, t) = -nh (x, t) [F(x, t) - B(x, t)] , 

where F and B are scalar functions representing the for- 
ward and backward fluxes in the PhC. Substituting in 
Eq. © gives 



„ , „ n(x, t) „ n* „ 
F x {x,t) + ^^F t + ^F 
c c 

B x (x,t)-^B t -^B 
c c 



_ 1 fnt _ nx 
2 V c n 
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B. 
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We rewrite the system in matrix form as 
nW t = MW, 



(5) 



where 



F 



w=[ B , M = 



-cd x - n t 



(^ + lf) ^ x -n t 



Near a gap, the solution of Eq. <j5j) consists predominantly 
of two spectral components on each of the two bands. We 
write each component as the product of a carrier uni- 
directional Floquet-Bloch (FB) mode ^ f/b an d a slowly 
varying envelope f /b, respectively, i.e., 



W = [f( Xl , h)* f (x ) + 6(xi, ti)*6(!ro)] e" 



(6) 



where lu c = 2irc/\ c and 2//6 arc the eigenmodes of 
Eq. ([5]). We also formally distinguish between fast varia- 
tions of the carrier waves (through xq and to) and slower 
variations of the envelopes (through x\ and t\). 

Assuming the central frequency of the input pulse is 
tuned exactly to the crossing point, substituting Eq. flS]) 
in Eq. multiplying by and Vf^, respectively, and 
integrating over xq allows us to remove the dependence 
on the fast scales. Then, for a piecewise-uniform mod- 
ulation, we obtain the following equations for the time 
evolution of the forward/backward envelopes 



1 + Am^'mff) J f t (x, t) + v g f x + Am^ (m t - iw c m(t)) / 
= Am {od) (iu c m(t)b - m t b - m(t)b t ) , 



(7) 



and 



(l + Am (d) m(i)) b t (x, t) - v g b x + Am (d) (m t - iw c m(t)) b 
= Am^ (iw c ro(*)/ - m t f - m(t)/ t ) , 

where v g = cj ' ■ s Jn\ri2 is the group velocity and mS d > 

77j/ b * with 



(8) 



m ff = m b \ m (od) _ m fb 



^}{x)p(x)^_Jx)dx. 



For simplicity of notation, we have dropped the subscript 
of the slow scales coordinates. 

In Eqs. (O-©, the terms proportional to m(t) on the 
left-hand-sides describe the correction to the wave veloc- 
ity, or cquivalcntly, the correction to the phase accumula- 
tion, induced by the modulation. The terms proportional 
to rrit(t) represent the change of the electromagnetic en- 
ergy induced by the modulation. The terms on the right- 
hand-sides of Eqs. (|7|)-(|8|l are responsible for reversal cou- 
pling, i.e., to the vertical transitions lu + 5lj — > u> — 5u>, 
see Fig. lb. One can interpret the magnitude of the 
associated coefficients m(° d ) as the measure of the phase 
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mismatch between the modes immediately above and be- 
low the crossing frequency. The coupling tends to zero 
as the index contrast approaches zero, as expected. In 
all other cases, the dependence of the coupling coefficient 
on the indices is non-trivial. 

Eqs. Q-JH]) can be solved analytically if we neglect 
all the coupling terms in Eq. (JT]). This is justified as 
long as the backward component, which is absent prior 
to the modulation, is small compared with the forward 
component. Under this assumption of a weak coupling, 
the solution of Eq. ([7]) can be shown to be 



Am (d) 



iuj c m{t') — mt(t') 
l + AmWm(f) ' 



where x^ 



dt' 



represents a frame 



1+Am( d )m(i') 

moving with the forward pulse and fo is the envelope of 
the pulse traveling in the forward direction before the 
onset of the modulation. In a similar manner, it can be 
shown that the backward wave ([8]) is given by 

= Am^e^%{x^,t), (9) 



where 



5 (*«*)= J h(t')f (>) - J - 



2v g dt» , 



and 



hit) 



+ Amt^mft") 
iuj c m(t) + mt{t) 



[l + AmWm(i)]' 



dt' 



is a frame moving with 



where x^ — a 

the backward pulse; this shows that the wave-front has 
indeed been reversed. At times long after the modulation 
has ended, b is effectively given by a convolution of the 
forward wave fo with h(t), which thus plays the role of 
the impulse response of the system. Since the durations 
of fo and h(t) are T mo d and T p , respectively, Eq. © 
shows that accurate reversal can be achieved only in the 
non-adiabatic limit 



(10) 



in which case the impulse function is essentially a delta 
function. In order to see that more clearly, consider a 
Gaussian modulation and Gaussian pulse, 



m(t) 



fo(x-v g t) = e °i T i . (11) 



Then, for a small modulation (A <C 1), we get that 
x (f> b ) = x zp Vg f^ anc j the backward wave amplitude be- 
comes 

b = -iMTe/zAmMe^O/o ( ^±^) , (12) 



T, 



eff 



with T eff = T mod T p / ^AT? nod 



7^2 

J-p . 



Thus, the reversal 

is accurate if condition (|10[) is satisfied; otherwise, the 
reversed pulse is broadener than the input pulse. 

Eq. (O can also be solved asymptotically without any 
simplifying assumptions. However, whenever the back- 
ward wave becomes comparable in magnitude to the for- 
ward wave, the weak coupling solutions (|9|) or (|12|) are 
not valid anymore. In such cases, Eqs. ([Z])-® should no 
solved without approximation. 

We now turn to discuss some practical aspects of a 
Reversal Mirror (RM) design for electromagnetic pulses. 
First, the RM should be long enough to contain the pulse 
during the modulation. Simple estimates show that a few 
cm long RM can contain and reverse pulses not longer 
than a few tens of picoseconds. Second, the thicknesses 
of the layers determine the carrier (crossing) frequency 
according to Eq. ([2]). Third, with regards to modula- 
tion techniques, picosecond and longer modulation times 
can be achieved by exploiting the electro-optical effect 
or using carrier injection in semiconductors. In the lat- 
ter case, index modulations of the order of 1% at about 
100GHz [25[ or even at lps-scale [28| were achieved in 
silicon devices. Shorter modulations must rely on an in- 
tense pulse in a pump-probe configuration, see e.g., fl9| . 
In these cases, however, the reversal efficiency will be 
smaller due to the shorter modulation time, see Eq. (|12l) . 
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FIG. 3: (Color online) (a) Amplitude of an asymmetric pulse 
(Eq. (21; red line) and the associated backward wave envelope 
(Eq. J5J; blue line) at the input side of the QWS as a function 
of time, (b) A spatio-temporal countour map of forward and 
backward envelopes in (a). 

Finally, we demonstrate the performance of our RM 
through numerical simulations. In Fig. E)Ja), we plot the 
wave amplitude at the input side of the PhC as a func- 
tion of time. The reversed pulse has a somewhat lower- 
amplitude, but the leading and trailing edges have clearly 
exchanged roles. We also show that the solutions of the 
wave equation @ and the envelope equations (O-© arc 
in excellent agreement. Fig. [3th) shows a spatio-temporal 
countour map of the pulse propagation. It shows that the 
modulation, occuring once the pulse is in the middle of 
the RM, causes a fast and complicated dynamics after 
which the pulse splits into a reversed and a (somewhat 
delayed) forward component. 

A comparison of the reversal efficiencies as a function 
of the modulation strength, which also includes the an- 
alytical approximation (|12[) , is shown in Fig. @] We em- 
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FIG. 4: (Color online) (a) Reversal of a T p = lOOfs, unit- 
amplitude Gaussian input pulse under a Gaussian modula- 
tion as a function of the index change A for ni = 3.45, 
712 = 1, Ac = 1550nm, and T moc i = lOfs. Shown are numerical 
solutions of the wave equation (Eq. (|4|; blue dots) vs. the 
solution of the envelope equations (Eq. flBjl: black circles) and 
the analytical solution (Eq. {12}; red solid line), (b) Same as 
(a) for T mo d — lps and T p = lOps. 

ploy realistic parameters corresponding to a silicon-air 
QWS with A c = 1550nm; such a PhC can be fabricated 
on-chip with current technology pol ]. First, we perform 
simulations of ~ 20-cycle pulses. Fig. HJa) verifies that 
all three solutions are in good agreement, thus, validat- 
ing the analysis. Fig. HJb) shows simulations for longer 
pulses for which memory and running time required for 
the solution of Eq. Q are very long. Thus, we only 
show the solutions of the envelope equations (0-© vs - 
the analytical solution. As in Fig. Hta), there is very 
good agreement between the numerical and analytical 
solutions up to high efficiencies (~ 50%). At even higher 
efficiencies, the analysis overestimates the reversal effi- 
ciency given by the solution of the envelope equations. 
This is because at such high efficiencies, the forward wave 
amplitude is significantly decreased, so the forward-to- 
backward wave coupling decreases as well. Nevertheless, 
the solutions of the envelope equations show that a 100% 



reversal efficiency can be achieved with index modula- 
tions only slightly stronger than those predicted analyt- 
ically. Thus, our RM has comparable performance to 
the previously suggested time-reversal schemes in index- 
modulated coupled-resonator arrays of optical waveg- 
uides [U [22| and wave-mixing- based schemes [l9| . How- 
ever, in contrast to these schemes, the QWS can admit 
and reverse pulses of very broad spectrum. In addition, 
our system does not suffer from a deterioration of perfor- 
mance due to non-adiabatic modulations. 

In summary, we have shown how to time-reverse wave 
packets by exploiting the unusual band structure of zero- 
gap PhCs using index modulations that can be very weak 
and much slower than the wave period. The suggested 
design was shown to allow for a 100% reversal efficiency, 
and can be applied to pulses of unprecedented broad 
spectrum. Our design can be fabricated with contem- 
porary technology, e.g., on a silicon chip. Our RM docs 
not require any prior knowledge on the pulse properties 
and can be implemented in almost any spectral range 
and pulse duration. Higher-dimensional realizations of 
zero-gap PhCs, e.g., having a hexagonal (graphene-like) 
symmetry, may open the way to additional reversal func- 
tionalities and may even be used as a building block in 
more sophisticated time-reversal schemes or other wave 
manipulation schemes like light stopping, slowing-down 
etc.. Finally, since our scheme does not rely on any con- 
cept which is unique to optics, it can be employed in 
other zero-gap systems for which time-reversal may not 
have been accessible so far, such as quantum systems [l7| . 
spintronics Q and microwave transmission- lines [f| . 

Y.S. would like to thank A.I. Fernandcz-Domingucz 
and A. Aubry for innumerable valuable discussions. Y.S. 
was supported by the Royal Society International New- 
ton Fellowship. J.B.P was supported by the EU project 
PHOME (Contract No. 213390). 



[14] 



J. Joannopoulos et al., Photonic Crystals (Princeton Uni- [15 

versity Press, 2008), 2nd ed.. [16 

R. Olsson and I. El-Kady, Meas. Sci. Technol. 20, 012002 [17 

(2009). [18 
A. Geim, Science 324, 1530 (2009). V.V. Cheianov et al. 

Science 315, 1252 (2007). [19 

J. Pendry, Science 306, 1353 (2004). [20 

M. Antoniades and G. Eleftheriades, IEEE Antennas and [21 

wireless propagation letters 2, 103 (2003). [22 

M. Berry, P. Roy. Soc. A-Math. Phy. 462, 1629 (2006). [23 

E. Rashba, J. Supercond. 15, 13 (2002). [24 
O. Peleg et al, Phys. Rev. Lett. 98, 103901 (2007). 

M. Fink, J. Phys. D 26, 1333 (1993). [25 

D. Pepper, Laser handbook, Vol. 4, p. 333 (North- [26 
Holland Physics, Amsterdam, 1988). 

G. Lerosey et al, Science 315, 1120 (2007). [27 

G. Lerosey et al, Phys. Rev. Lett. 92, 193904 (2004). 

J. Pendry, Science 71, 322 (2008). A. Aubry and J. [28 

Pendry, J. Opt. Soc. Am. B 27, 72 (2010). 

X. Li and M. Stockman, Phys. Rev. B 77, 195109 (2008). 



Z. Yaqoob et al, Nat. Phot. 2, 110 (2008). 

A. Ruffin et al, IEEE J. of Qu. Electron. 38, 1110 (2002). 

F. Cucchietti, J. Opt. Soc. Am. B 27, A30 (2010). 

A. Yariv, D. Fekete and D. Pepper, Opt. Lett. 4, 52 

(1979). 

D. Miller, Opt. Lett. 5, 300 (1980). 

D. Marom et al, Opt. Lett. 25, 132 (2000). 

M. Yanik and S. Fan, Phys. Rev. Lett. 93, 173903 (2004). 

S. Longhi, Phys. Rev. E 75 026606 (2007). 

O. Kuzucu et al, Opt. Exp. 17, 20605 (2009). 

M. Notomi, Rep. Prog. Phys. 73, 096501 (2010). S. Fan 

et al, Opt. Phot. News 18, 41 (2007). 

M. Lipson, J. Lightwave Technol. 23, 4222 (2005). 

P. Yeh, Optical Waves in Layered Media (Wiley- 

Interscience, 2005), 2nd ed.. 

C. de Sterke, D. Salinas and J. Sipe, Phys. Rev. E 54, 
1969 (1996). 

A. Chin et al, Appl. Phys. Lett. 69, 653 (1996). 



